Dynamical Variational Principles for Strongly 
Correlated Electron Systems 



Michael Potthoff 

Institut fur Theoretische Physik und Astrophysik, Universitat Wurzburg, Am 
Hubland, 97074 Wurzburg, Germany 

Abstract. The self-energy-functional approach (SFA) is discussed in the context 
of different variational principles for strongly correlated electron systems. Formal 
analogies between static and dynamical variational approaches, different types of 
approximation strategies and the relations to density-functional and dynamical 
mean-field theory are emphasized. The discussion elucidates the strengths of the 
SFA in the construction of new non-perturbative approximations but also the lim- 
itations of the approach and thereby opens up future perspectives. 

Variational approaches have a long and successful tradition in the theory 
of condensed-matter systems as they offer a smart, controlled and system- 
atic way to treat the problem of electron-electron interaction. A well-known 
variational approach is Hartree-Fock (HF) theory. It is based on the Rayleigh- 
Ritz principle and provides a practicable and consistent mean-field descrip- 
tion of an interacting electron system. As quantum fluctuations are neglected 
completely, HF theory must be classified as a static mean-field approxima- 
tion. This may be contrasted with dynamical mean-field theory (DMFT) jll2j 
which includes temporal fluctuations in the mean-field picture. The DMFT, 
however, cannot be derived from the Ritz principle. It must be constructed 
from some dynamical variational principle which involves a dynamical (i.e. 
time- or frequency-dependent) quantity as the basic object. Dynamical varia- 
tional principles have already been suggested in the sixties |3I4| but, compared 
to the Ritz principle, were employed with rather limited success only. This 
brings up the following questions: What are the similarities and the differ- 
ences between different variational principles and approximation strategies? 
How can the DMFT be considered as an approximation within a variational 
concept? Can dynamical variational principles be used for constructing prac- 
ticable and non-perturbative approximations different from the DMFT? An 
attempt to answer these questions straightforwardly leads to the self-energy- 
functional approach (SFA) 5] suggested recently. The purpose of this paper 
is to discuss different variational approaches and to place the SFA into this 
context with the objective to explore possible future developments. 

1 Variational Principles and Approximation Strategies 



Consider a many-electron system in the volume V, at temperature T and 
with chemical potential fi. It is characterized by a Hamiltonian Ht.u = 
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Ho(t) + Hi(U) consisting of a one-particle and an interaction term Hq and 
Hi , respectively, which depend on the "model parameters" t and U (a matrix 
notation is used). In second-quantized form, 



where an index (e.g. a) refers to the states of a one-particle basis. 

The characteristic of a variational approach is a certain physical quantity 
X to be varied, as e.g. the statistical operator, the electron density, the (local) 
Green's function, the self-energy etc. Clearly, at equilibrium this quantity will 
depend on the model parameters: Xt } u (and on V,T and /i). 

In a variational approach, the quantity is considered to be a variable. The 
first task is to express a thermodynamical potential fl (the grand potential, 
for example) as a function(al) of this variable: Qt u\X\. As this functional is 
characteristic for the model system {Q , it will depend on the model parame- 
ters. At the equilibrium or "physical" value, i.e. at X = Xt,u, one must have 



Sh,v[Xt,u\ = fk,u where f2 t<u = -Tlntrexp(— (iJ (*) + H 1 (U) - fj,N)/T). 



Furthermore, the functional Cl t ,u[X\ should be constructed such that 
it becomes stationary at the physical value: 8f2t,u[X = ^t,u] = 0. This 
variational principle determines Xt t u once the functional is known. Note that 
the domain of the functional must be specified in addition since in most cases 
Xt,u must satisfy some constraint or sum rule or normalization condition. 

Even if the functional is known, however, it is usually impossible to eval- 
uate f2t,u[X] for a given X, and one has to resort to approximations. One 
may distinguish between three types of approximation strategies: 

In a type-I approximation one derives the Euler equation 6f2t,u[X]/6X = 
first and then chooses (a physically motivated) simplification of the equa- 
tion afterwards to render the determination of Xt : u possible. This is the 
most general but worst type, as normally the approximated Euler equation 
no longer derives from some approximate functional. This may result in ther- 
modynamical inconsistencies. 

A type-II approximation modifies the form of the functional dependence, 
^t,u[' • •] — * &t,u[' • •], to get a simpler one that allows for a solution of the 
resulting Euler equation 6f2t,u[X]/6X = 0. This type is more particular but 
yields a thermodynamical potential consistent with Xt.u- Generally, how- 
ever, it is not easy to find a sensible approximation of a functional form. 

Finally, in a type-Ill approximation one restricts the domain of the func- 
tional. The restriction comes in addition to those conditions that are physi- 
cally necessary anyway (e.g. normalizations) and requires a precise definition 
of the domain. This type is most specific and, from a conceptual point of 
view, should be preferred as compared to type-I or type-II approximations as 
the exact functional form is retained. In addition to conceptual clarity and 
thermodynamical consistency, type-Ill approximations arc truely systematic 
since improvements can be obtained by an according extension of the domain. 
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Note that any type-Ill approximation can always be understood as a type-II 
one (and type-II approximations as type-I) but not vice versa. 

2 Various Variational Approaches 

In the following different variational principles and possible approximations 
shall be discussed according to this scheme, starting with Hartree-Fock and 
density-functional theory as illuminating examples. 

2.1 Ritz Variational Approach 

In the Ritz variational approach the ground-state energy is considered as a 
functional of the quantum state |<P). There is a generalization of the Ritz 
principle to arbitrary temperatures by Mermin [S] . Here the basic variable is 
the statistical operator p characterizing the system's (mixed) state, and the 
grand potential as a functional of p reads: 



Following the classical calculation of Gibbs, it can easily be shown |H] that 
on the set of normalized and positive definite operators, stationarity of the 
functional © is achieved for p = p t .u = e -{H t ,u-nN)/7 ^^ e -(H t ,u-^N)/Ty 
One also has f2t.u[pt,u] — ^t,u- An additional feature of the functional J5J 
consists in the fact that Qt,u[p] > ^t,u[pt,u] for & n y P- This "upper-bound 
property" is extremely helpful but specific to the Ritz principle. 

For a many-electron system and an arbitrary p, the computation of the 
trace in Eq. (J2J is an exponentially hard problem. A nice type-III approxi- 
mation is the HF approach: Here the variational search is restricted to the 
subclass of disentangled statistical operators, i.e. statistical operators corre- 
sponding to independent-electron states. This can be made precise by intro- 
ducing the important concept of a reference system: 

A reference system is a system with a different (microscopic) Hamiltonian 
Hf u' 7^ *) U' 7^ 17) but with a macroscopic state characterized by the 
same values of the thermodynamic state variables as the original system : 
V = V, T' = T and p! = p. The sole purpose of the reference system is 
to specify the domain of the functional J2J): Trial statistical operators are 
taken from the reference system, p = pt',u'i an( i are varied by varying the 
parameters t' and U within a certain subspace. Hence, the choice of the 
reference system (and the parameter subspace) defines the approximation. 

The HF approximation is given by the choice U = and t' arbitrary, i.e. 
by trial states p t >.o = e _(H '' 0_AlAr) / T /Z t /. . Inserting into Eq. © yields 



Ot,u[p\ = te{p(Ht,u - pN + Tlnp)) . 



(2) 



fh,uW,o] = flt',o + tr(^ >0 (77o(t) + Hi(U) - H (t')) . 



(3) 



The remaining trace can be computed easily using Wick's theorem as pt'.o 
derives from a non-interacting Hamiltonian. The variational parameters t' are 
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fixed by the conditions 8f2t.u[pt' ,o]/8t' = 0. These are exactly the well-known 
HF equations as can be seen by some straightforward manipulations. 

One learns that type-Ill approximations can be constructed conveniently 
by the concept of a reference system. On the one hand, the reference system 
should comprise a large space of parameters t' and U to generate a powerful 
approximation. On the other hand, the parameter space must be restricted 
strongly to keep the calculations feasible. 

2.2 Density- Functional Approach 

For a many-electron system the statistical operator or, at T — 0, the ground- 
state wave function actually is an object that is by far too complex. The 
relevant information is much more efficiently stored in integral quantities, 
such as the electron density. This is the starting point of density-functional 
theory (DFT) |7I8I9| . The density, i.e. the quantum-statistical average of the 
one-particle density operator n(r) = tr(pn(r)), represents the basic variable. 
Normally DFT aims at the inhomogeneous electron gas at T = but can 
also be applied to discrete lattice models TO] and finite temperatures UJ. 

The grand potential f?t,u obviously depends on the model parameters. 
Due to the Hohenbcrg-Kohn theorem [7], however, it can also be considered 
as a functional of the density n which is stationary at the physical density: 
8Qt.u[n] = for n = nt.u- Furthermore, if evaluated at n = nt y u, it 
yields the exact grand potential: fit,u[nt.u] = ^t,u- Keeping the notations 
introduced above, n is a matrix with n a p = tr(p c^c^), and the functional 
reads (cf. Refs. HO]]): 



Here the trace refers to the one-particle orbitals a, and .F[/[n] is a universal 
functional, i.e. it depends on the interaction parameters only. Using the Kohn- 
Sham idea |8I1()| , the resulting Euler equation has the form of a one-particle 
Schrodinger equation. 

The variational principle 8fit,u[n\ = is rigorous but cannot be eval- 
uated as Fu[n] is generally unknown (after separating the Hartree and a 
kinetic-energy term, the remaining exchange-correlation functional is not ex- 
plicit). Due to the universality of Fj/fn], however, the density- functional for 
a reference system with modified one-particle parameters t' reads fit',u[ n ] = 
tr(i'n) + and thus f2t,u[ n t' ,u] — ^t',u + tr((i — t') n t',u") which can 

be exploited for a type-Ill approximation. Choosing as a reference system 
Hf ,u the homogeneous electron gas, however, turns out to be too restrictive, 
as this implies a spatially constant density. The local density approximation 
(LDA) |7l8j . on the other hand, has proven to be very successful 9 . At least 
for systems with weakly varying density it is well justified. The LDA, how- 
ever, is no longer a type-Ill approximation but a type-II one as the form of 
the (exchange-correlation part of the) functional Fjj [n] is approximated to 
have a local dependence on the density only. 




(4) 
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As the proof of the Hohenberg-Kohn theorem is based on the Ritz prin- 
ciple [7j , the upper-bound property is transferred to the exact functional Q , 
i.e. fit : u[n] > ^t,u for any n, but is lost within the LDA due to the type-II 
character of the approximation. 

2.3 Time-Dependent DFT 

The weak point of the DFT consists in its inability to describe excitations 
(see, however, Ref. QU). This is due to the fact that the Hohenberg-Kohn 
variational principle is built on the static electron density. Information on ex- 
citation properties is contained in dynamic response functions which are ac- 
cessible in principle via time-dependent density-functional theory (TD-DFT) 
|12| . In TD-DFT one considers a situation with a time-dependent Hamilto- 
nian and focuses on the time-dependent density n(r,t) = (!f - (t)|r';,(r)|S'(t)) 
resulting from a solution \^(t)) of the time-dependent Schrodinger equation 
as the basic variable. Here the action A = J dt(^(t)\id/dt — H(t)\&(t)) can 
be understood as a functional of n(r, t), 



where t(t) are explicitly time-dependent one-particle parameters. Contrary 
to usual DFT, the variational principle SAt t u[ n ] = does not derive from the 
Ritz principle, and consequently there is no upper-bound property available. 
Type-II approximations can be constructed by approximating the universal 
but unknown part Bjj[n] of the functional to make it explicit. Far from 
equilibrium, however, there is no general recipe. 

2.4 Dynamical Variational Principle 

In the linear-response regime close to equilibrium, excitations are described 
by Green's functions. The one-electron Green's function G a p(oj) — ((c a ; Cg))^ 
is the basic quantity in the dynamical variational approach of Luttinger, 
Ward, Baym and Kadanoff |3I4| . Employing a coupling-constant integration 
H3, the grand potential can be understood as a functional of G: 



where Tr = T£„e M "° tr and u>„ = (2n + 1)ttT are fermion Matsub- 
ara frequencies. Furthermore, Gt,o = {oj + /i — t) -1 is the U = (free) 
Green's function and <Pu [G] the (universal) Luttinger- Ward (LW) functional 
defined as the sum of all dressed closed skeleton diagrams [Hj- By construc- 
tion, f2t.u[Gt.u] = ^t.u- In arbitrary order in perturbation theory one has 
6$u[G]/6G = TSu[G}. Therewith, the Euler equation 8fl t ,u[G]/5G = is 
given by G _1 — G7q + £u[G] = which is Dyson's equation. This shows 
that flt,u[G] is stationary at the physical Green's function G = Gt,u- 




to 



(5) 




(6) 
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The LW functional is formally given by a diagrammatic sum that cannot 
be carried out in practice. A self-evident type-II strategy is to sum up a 
suitable subclass of diagrams to obtain an approximate but explicit expression 
for <Pjj [G] . In this way the HF approximation can be recovered but there are 
also new approximations like the fluctuation-exchange approximation [13114) . 
These "conserving approximations" , however, are necessarily restricted to the 
weak-coupling regime. 

A type-Ill approximation, on the other hand, would be non-perturbative 
by construction. Consider a reference system with modified one-particle pa- 
rameters: Hf.u = Ho(t') + H\{U), This defines the domain of the func- 
tional © to consist of Green's functions Gf,u with arbitrary t' . The in- 
teraction is kept fixed (U' = U). To evaluate the functional (JBJ) at Gt>,u 
requires the evaluation of ^u[Gf,u], m particular. Due to the universal- 
ity of $u[- • •] (no t dependence) and due to the choice U' = U ', one has 
$u[G t >,u] = fi VtU - TrlnG^.t/ + Tr{(G^ - G^G^u). Thus, 

fitMGt',u] = Shr,u - Ti(G^G v>u ) + Tr(G^ a G ViU ) . (7) 

Hence, on any domain specified by a suitable subspace of one-particle pa- 
rameters t' which renders the solution of the reference system possible (for 
fixed U), the functional JHJ) can be evaluated exactly. A possible (but over- 
simplified) example is the choice t' = 0. It reduces the reference model to the 
atomic limit where the computation of Gt> ,u an d fit' .u m is easy. Cluster 
approximations represent straightforward generalizations of this example. 

Unfortunately, this type-Ill approach for Eq. © yields nothing new: Since 
Tt{G^o - Gt\)G t >u = tr(< - t')n v>v with the one-electron density of the 

reference system n v >v = TJ2 n e lu) "- 0+ G t > ,u{iu n ), one gets Qt,u[Gt',u] = 
^t,u[pt' ,u] with f2t,u[p\ given by Eq. (j2J), i.e. the same as in the Ritz varia- 
tional approach. Interestingly, this implies that upper bounds for the grand 
potential can be obtained, i.e. fit,u[Gf ,u] > ^t,u for arbitrary t'. 

2.5 Dynamical Mean-Field Approach 

Equipped with these insights, one can address the question of deriving the 
DMFT from a variational principle. Originally, the DMFT was introduced as 
the exact theory of lattice models with local (Hubbard-type) interactions in 
infinite spatial dimensions D = oo . Later on, it was recognized [15116) that 
in D = oo the lattice model Ht,u can be self-consistently mapped onto an 
impurity model Hf .u with the same interaction U . Using this self-consistent 
mapping procedure as an approximation ( "dynamical mean-field approxima- 
tion"), one can treat lattice models for any finite D. 

Instead of considering Dyson's equation in the form, G = (G^q — [G] ) ~ 
(with Su[G}) = (1/T) S$u[G])/SG), which is solved by the exact G t<U , the 
DMFT considers the following simplified equation between local quantities at 
lattice site i: 

(G) u = {G- t l - Eu[Q\)« l . (8) 
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Here Su[G] is the derivative of the LW functional but with local vertices 
only as it is the case for an impurity model. Clearly, this is a type-I ap- 
proximation. Eq. (JHJ is often called the DMFT self-consistency condition. 
This is because its solution is achieved by an iterative procedure in prac- 
tice: Starting with a guess for S, one computes the local lattice Green's 
function as (G)u = (G^q — -S 1 )^ 1 at first. This is not yet a solution of Eq. 

JSJl since in general £ ^ Sjj[G] for this G. For the necessary update of 
S define (Gf,o)u = 1/(1/(6?)" + (S)u). Assuming that Gf,o can be un- 
derstood as the free impurity Green's function of an impurity model Ht'.u 
for some t' , the (numerical) solution of the impurity problem yields a new 
£ = £xj[(G)u] — Sjj [G]. Iteration of this cycle until self-consistency yields 
a solution G of Eq. JHJ. Note that the resulting DMFT self-energy is local. 

Here, the question is whether Eq. JHJ can be interpreted as an Euler 
equation of some variational principle. Starting with the functional ©, one 
can try a type-II approximation by replacing <Pjj [G] with the LW functional 
of the impurity model <% [G] . This implies (l/T)S$u[G]/6G = S V [G\, and 
the resulting Euler equation reads: G = G^q — Su[G]. This equation 

is easily seen to be equivalent with Eq. JSJ since £u[(G)u] = Sjj[G\ by 
definition. Hence, DMFT can be understood as a type-II approximation. 
Another functional has been suggested recently |17|: 

Q t ,u[G] = Trln 1 - Tr(Z u[G]G) + <t>u[G] . (9) 

^tfi ~ ^u[G\ 

Clearly, ftt.u[Gt,u] = &t.u, and furthermore the corresponding Euler equa- 
tion, ((G^o - SuiG})' 1 - G)-(6Su[G]/SG) = 0, is equivalent with Dyson's 
equation, G = (G^q — SuiG])^ 1 (assuming local invertibility of the func- 
tional Su [G]). The functional therefore yields a valid variational prin- 
ciple. As a type-II approximation, one may replace ^c/[G] — + <Pu[G] and 
Su[G] -> Su[G] = (l/T)S$u[G]/5G in the functional ©. The resulting 
Euler equation is equivalent with the DMFT self-consistency equation JHJ 
which implies that DMFT can also be understood as a type-II approxima- 
tion to the functional ©. 

Attempts to prove that a stationary point of the type-II approximated 
functionals © or JHJ must be an extremum have failed |17| . Furthermore, 
while (as shown above) a type-Ill approximation to the principle based on 
Eq. © with U fixed yields upper bounds to the grand potential, the DMFT 
cannot be obtained as a type-Ill approximation starting from Eqs. © or @: 
Choosing an impurity model Ht\u as a reference system to generate trial 
Green's functions and to define a restricted domain of the functional © or 
jnjl, respectively, concurrently means that the optimal Green's function will 
be local. This is obviously a very poor approximation for the Green's function 
of a lattice model and differs from the DMFT result. The discussion shows 
that the question whether or not the DMFT grand potential is an upper 
bound to the true grand potential is still open. 
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2.6 Functionals of the Local Green's Function 

It is also possible |18I19| to focus on the local Green's function G*- 100 -* = (G)u 
(instead of the full G) and to set up a variational principle of the form 

sn tl u[G (loc) ] = o- (io) 

A functional which is stationary at the physical G^ 00 -* = G^^j and which 

yields Qt,u\G^ ^ ] = ^t,u can be constructed order by order in the interac- 
tion strength |18| . Unfortunately, the diagrammatic formalism is much more 
cumbersome as compared to the construction of the LW functional. As is 
shown in Ref. |19| . the dynamical mean-field approximation is equivalent 
with a simple ([/-independent) approximation to the kinetic-energy part of 
the functional. So the DMFT appears as a type-II approximation again. 



2.7 Self-Energy-Functional Approach 

The motivation to characterize the DMFT as a type-Ill approximation is 
the following: If it is possible to recover the DMFT merely by restricting 
the domain of the functional corresponding to an exact variational principle, 
different choices of the domain will place the DMFT in a systematic series 
of different and possibly new approximations which, as the DMFT, are all 
non-perturbative and thermodynamically consistent. 

For this purpose it is helpful to focus on the self-energy. Within the self- 
energy- functional approach (SFA) [5], the self-energy functional 

n tiU [S]=Trla 1 +Fu[Z] (11) 

is considered. Here, F V [E] = <Pu[G[£}} - Tr (S Gu[S}) is the Legendre 
transform of the Luttinger-Ward functional which is well defined provided 
that the functional Su[G] is invertible locally. is universal (indepen- 

dent of t) by construction and -(l/T)8Fu[S]/8E = G V [E] which is the 
inverse of the functional £u[G]. Obviously, Slt,u[£t,u] = ^t,u- The Euler 
equation 8Q t ,u[^]/ 8S = is given by (G^q — S)^ 1 = G[S] and equivalent 
with Dyson's equation. 

To construct a type-Ill approximation, a reference system Hf.u = Ho(t')+ 
Hi(U) with unchanged interaction part is considered. The one-particle pa- 
rameters t' are taken such that the different "correlated" sites (non-zero on- 
site interaction) are decoupled. Instead, t' shall include an arbitrary hopping 
to "bath" sites (zero on-site interaction) with arbitrary one-particle energies. 
In case of the Hubbard model on a lattice with L sites, the corresponding ref- 
erence system constructed in this way is a set of L decoupled single-impurity 
Anderson models (which in case of translational symmetry are equivalent). 
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Trial self-energies St',u are local by construction. The Euler equation result- 
ing from this type-Ill approach reads df2t,u[^t',u]/dt' = 0, i.e.: 

((G^ 1 - S t ,.u)- 1 G v [E V ,u\) ■ d -^f L = . (12) 

Now let t' (the bath parameters) be such that Gf : u solves the DMFT self- 
consistency condition JSJ. Since Uu[Gf,u] — ^t',U, ° ne bas {Gt>.u)u = 
(G^Q - Et^u)^ 1 . Hence, this t' solves Eq. JT2}. (Note that dE t ',u/dt' is 
local). So by a restriction of the domain of the self-energy functional to 
local self-energies, the DMFT is characterized as a type-Ill approximation. 

Interestingly, a type-II approximation does not work: A replacement of 
the form Fu[S] — » Fu[S] in Eq. ifTTl) yields the Euler equation Gu[S] = 
(G^ 1 - i;)" 1 where Gu[S] = -{l/T)SFu[E]/5E. If this was equivalent 
witb the DMFT self-consistency condition, a local self-energy would be a 
solution. This would imply, however, that Gu[S] is non-local for a local S. 
Hence, Fu [S] cannot be the Legendre transform of <Pu [G] where <!>u [G] (as 
above) is the Luttinger-Ward functional with vertices restricted to a single 
site. An alternative choice for F[/[27] ; however, does not suggest itself. 

One can conclude that a functional of the self-energy is necessary and 
sufficient to obtain the DMFT as a type-Ill approximation while a functional 
of the Green's function is necessary and sufficient to get the DMFT as a 
type-II approximation. The decisive point is that rather a local self-energy 
can be tolerated as an approximation than a local Green's function. 

3 New Approximations 

The immediate return of these considerations is a number of non-perturbative 
and thermodynamically consistent type-Ill approximations as shown in Fig. 
^ These differ from the DMFT by a different restriction of the domain for the 
self-energy functional (|llf> . i.e. by a different reference system with a differ- 
ent subspace of variational parameters t' but the same interaction (U = U). 
The evaluation of a type-Ill approximation requires the repeated computa- 
tion of the grand potential and the Green's function or self-energy of the 
reference system to get f2t,u[Et>,u] which must be optimized with respect 
to t' subsequently. 

The class of possible reference systems is essentially spanned by two pa- 
rameters, namely n s — 1 which is the number of additional bath sites per 
correlated site and N c which is the number of correlated sites in a cluster 
that is decoupled from the rest of the correlated sites (Fig. The DMFT 
is obtained for N c = 1 and n s = oo since a continuous bath (n s = oo) is 
necessary to represent an arbitrary local free Green's function. 

The choice N c = 1 but n s < oo yields new approximations ("dynamical 
impurity approximations" , DIA) which are inferior as compared to the full 
DMFT but allow for much simpler and faster calculations when n s is small. 
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The most simple but non-trivial approximation (n s = 2-DIA) has been shown 
|5I2(J| to already cover the essence of the DMFT scenario [T] for the Mott 
metal-insulator transition in the Hubbard model. At the critical point for 
T = the calculations can be done even analytically ]20] > an d with increas- 
ing n s the grand potential, static quantities and the entire phase diagram 
rapidly converge to the full DMFT results |5I21| . The DIA is similar but su- 
perior as compared to the exact-diagonalization approach . Even for small 
n s the approach is thermodynamically consistent and, off half-filling, respects 
the Luttinger sum rule, for example. The DIA has also been employed suc- 
cessfully to study the influence of phonons on metal-insulator transitions in 
the Holstein-Hubbard model |22I23| . 

Nothing new is obtained for n s — oo and N c > 1: Here the SEA recovers 
the cellular DMFT [21]. (Note that the dynamical cluster approximation 25 
is a type-II approximation). More interesting is the case n a = 1 and N c > 1 
which turns out |2(i| to represent a variational generalization of the cluster- 
perturbation theory [27]. This V-CPT is well suited to describe phases with 
spontaneously broken symmetry and has been employed to study one-particle 
excitations and antiferromagnetic order in the D = 2 and D = 1 Hubbard 
model at half-filling [2H] and charge ordering in the extended Hubbard model 
|29| . A further application concerns antiferromagnetism in quarter-filled lad- 
der systems |3U]. An impressing example of the power of the V-CPT approach 
has been given recently in a study of the competition between antiferromag- 
netism and d-wave superconductivity in the hole- and electron-doped Hub- 
bard model [sH]- The question of phase separation is addressed in Ref. |H2j . 

Summing up, the SFA is able to unify different cluster theories and lo- 
cal approximations within a single and consistent framework which offers a 
large flexibility in the use of bath sites, ficticious fields, boundary conditions 
and particle reservoirs |26j . The formalism provides a controlled compromise 
between the demands for a non-perturbative and systematic theory working 
in the thermodynamic limit on the one hand, and the limited computational 
capabilities to diagonalize finite-size systems on the other. 
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4 Open Problems 

The self-energy-functional approach allows to construct a series of consistent 
approximations which improve systematically as N c — > oo. It is by no means 
clear, however, whether bath sites n s > 1 help to speed up the convergence 
with respect to N c and whether a cluster extension of DMFT or the V-CPT is 
more efficient. This can be clarified only empirically by considering different 
lattice models in different dimensions. As a few bath sites have turned out 
to be sufficient for reproducing the essential mean-field (N c = 1) physics in 
a number of studies of the single-band Hubbard model, further applications 
of the DIA are worthwhile to explore e.g. the mean-field phase diagrams of 
more complex (multi-orbital) models. Furthermore, one may also envisage the 
application of a simplified DMFT where a single (but continuous, n s = oo) 
bath is optimized for a multi-orbital model. This might be well justified for 
not too low temperatures. 

On the technical side, there are two main future tasks: The full diago- 
nalization and the Lanczos method which have been used so far, should be 
supplemented by a "reference system solver" based on stochastic techniques 
to improve the scaling of the numerical effort with the system size. Secondly, 
it would be advantageous to have an iterative technique at hand that directly 
yields a solution of the SFA Euler equation without the need for numerical 
differentiation. First results using full diagonalization |2J| are encouraging. 

On the conceptual side, the question for the possibility to give strict upper 
bounds to the grand potential is still open. Probably, a positive answer re- 
quires to establish a link to the Ritz variational principle. On the other hand, 
no example is known yet where the SFA grand potential at a stationary point 
is lower than the exact one. 

There are different directions into which the formalism may be extended. 
As the coherent-potential approximation for the disorder Anderson model has 
the same (mean-field) status as the DMFT for the Hubbard model, it sug- 
gests itself that a self-energy-functional approach can also be constructed for 
systems with disorder (and interaction). First applications |33| demonstrate 
that such a theory [31] is feasible. A challenge consists in the extension of the 
theory to include two-particle Green's functions in a generalized variational 
principle. Here the recently proposed functional-integral formulation of the 
SFA 35 can be helpful. Two-particle correlation functions are interesting by 
themselves and may furthermore facilitate an even greater flexibility in the 
choice of reference systems. At the same time such an approach could provide 
a conceptual clear way to treat models with non-local interactions. Currently, 
this problem is circumvented by a more pragmatic decoupling procedure |29|. 
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